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Presentation Outline 


justification - why write yet another turbomachinery code? 


approach - what does an object-oriented turbomachinery code look like? 


• results - how do I know the code works? 
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Justification 


there is still a need for 2-D design/analysis 
codes tend to be focused on one aspect 
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• specific, individual codes may have undesirable features 
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Problem Description and Assumptions 


CODE REQUIREMENTS: 

OTAC is applicable for 

- compressors and turbines 

- design and analysis 

- meanline and streamline 

- axial, centrifugal/radial, and mixed 


CODE ASSUMPTIONS: 

flow going through a blade row in an 
annulus from station (?) to station (2) : 

- steady-state, throughflow 

- circumferentially uniform 

- adiabatic, simple radial equilibrium 

- no change in mass flow rate 

- no streamline curvature 



ADDITIONAL GOALS: 

modular (loss models), good thermo, 
simulate unconventional architectures 
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OTAC Written in NPSS Environment 


• allows re-use of Numerical Propulsion System Simulation objects 

• model structure similar to NPSS engine cycle model 
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3-Stream OTAC Example Model 


Element object 



flow connection 



mechanical connection 
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FlowStation Object Extended from NPSS 


NPSS 1-D FlowStation (4 inputs): 

K P. 

MN 

rh 


OTAC FlowStation (7+1 inputs): 

K P. 

MN, a, <() 

rh 

radius 



+ relative frame angular speed: w 



National Aeronautics and Space Administration 


OTAC 


Streamtube in an Annulus 
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Multiple Streamtubes 
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NPSS BladeRow Objects 


the BladeRow represents the entire 
blade row and contains its own 
“sub-objects” 


each BladeSegment tracks a 
streamtube through a section of 
blade 
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entrance FlowStation 



exit FlowStation 


National Aeronautics and Space Administration 


OTAC 


9 




NPSS Solver 


Independents represent variables 
the NPSS solver is allowed to vary 


Dependents represent equations 
or conditions the NPSS solver must 
satisfy 
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FlowStation 

Independents 


< 


m 2 

ht2 

Pa 

a 2 


radius 2 


V MN 2 
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continuity 

^ml 



conservation of energy/Euler 

^t 2 ^ti = o)(r 2 V e2 -riV ei ) 

BladeRow j 


non-ideal process loss 

P t2 = P t2 ideal - AP t 

Dependents 


non-ideal process turning 

P2 — P blade "h ^ 



geometry constraint (radius) 

radius 2 — r maC f l i ne 



geometry constraint (area) 

Aflow2 ~ A machine ~ ^ blockages 
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Empirical Effects 


BladeRows contain Sockets, placeholders to insert code that calculates a 
certain variable such as non-dimensional pressure loss 



this allows for considerable versatility in applying losses to the simulation; 
other benefits include testing and proprietary considerations 
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Results 


comparison against other codes and calculations 


• investigation to determine even if the NPSS solver could reliably converge 
with matrix sizes over 50x50 


• more test cases have been run than shown here 
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Test Cases and Results 


• comparison of OTAC and HT0300 for a compressor IGV plus rotor, streamline, losses input 



Program HT0300, Richard M. Hearsey, 2011 
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Non-Dimensional Loss Parameter 


Test Cases and Results 


• comparison of OTAC and Ainley-Mathieson single stage turbine calculation, meanline, losses 
calculated 
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Rotor Incidence (degrees) 

A Method of Performance Estimation for Axial-Flow Turbines, D.G. Ainley and G.C.R. Mathieson, 1957 
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Adiabatic Efficiency 


Test Cases and Results 


• comparison of OTAC and HT0300 5-stage turbine calculation, streamline, losses calculated 
using Ainley-Mathieson with Kacker/Okapuu modifications 


Program HT0300, 

Richard M. Hearsey, 2011 


stage 4 results 
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Test Cases and Results 


• OTAC analysis of 5-stage turbine (from previous slide), streamline 
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Test Cases and Results 


• OTAC analysis of 2-stage compressor, streamline, losses calculated using Aungier correlations 



Corrected flow, Ibm/s 


Axial-Flow Compressors: A Strategy for Aerodynamic Design and Analysis, Ronald H. Aungier, 2003 
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Test Cases and Results 


• comparison of OTAC and Japikse & Baines centrifugal compressor calculation, meanline, losses 
input 


impeller exit 

OTAC 

Japikse 

Pt, psi 

31.17 

31.17 

Tt, R 

653.5 

653.7 

Vm, ft/s 

342.4 

342.4 

ve, ft/s 

843.8 

843.8 

P flow, degrees 

19.04 

- 19.04 

a flow, degrees 

67.91 

67.91 

slip factor 

0.8772 

0.8772 




diffuser exit 



Pt, psi 

30.04 

30.04 

Ps, psi 

26.71 

26.64 

a flow, deg 

55.99 

50.94 


Introduction to Turbomachinery, David Japikse and Nicholas C. Baines, 1994 
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Summary 


• OTAC proof of concept verified - correct results for compressors, turbines, 
axial, centrifugal, meanline, streamline, design and analysis 

• extensive work on turbine loss models: Ainley-Mathieson, Kacker-Okappu, 
Dunham-Came, Moustapha-Kacker-T remblay 

• compressor loss model based on Aungier’s method implemented 

• further work includes additional loss models, improved logic for choked flow 
operation 
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Backup Slides 
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Meanline BladeRow Equation Set 


continuity 

conservation of energy/Euler 
non-ideal process loss 
non-ideal process turning 
geometry constraint (radius) 
geometry constraint (area) 


rfl-m2 rfimi 

ht 2 - h t 1 = 0)(r 2 V e2 - r x l/ 01 ) 

P t2 = P t2 ideal - AP t 

P 2 — P blade ~i~ ^ 

syyf yyyf 

'2 ~ 'machine 

Afiow2 ^ machine ~ ^blockages 


note: at design, p blade and A machine may be 
input (direct-design) or varied to produce 
specific performance (indirect-design) 
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Streamline BladeRow Equation Set 


n 

n 

n 

n 

n-1 

1 


n-1 

1 


continuity 
energy/Euler 
loss condition 
flow follows blade 
geometry constraint 
geometry constraint 

1 dp Vn 

spanwise eq. - -j- = — 

p dr r 

geometry constraint 


n = number of streams 
i = stream number, 1 to n 
sum = aggregate value 


h t 2 i - h n i = ( 0 (r 2 .VQ 2 . - ri.Ve!.) 

^t2i P t2ideali 
P 2 i ~ P bladei 

sy* / 1/1 

2iiwi6r 2outeri 


sy* sy-1 

^ sultti iTLCLchine 


1 A pj _ v Q t 
Pi Ar i r t 

Afiow2 sum ~ A machine ~ A blockages 
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BladeSegment Object 


responsible for differences 
between certain flow states 
entrance 
exit - actual 
exit - ideal h, 
exit - ideal P t 
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BladeSegment Object 


multiple BladeSegments allow 
for radial variation of flow 
properties 
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BladeRow Object 


responsible for differences 
between BladeSegments 

holds blade row specific 
variables: annulus areas, 
number of blades, blade 
angles, power, etc. 
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Slide Master 
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